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SPB stars in the open SMC cluster NGC 371 
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ABSTRACT 

Pulsation in (3 Cep and SPB stars are driven by the k mechanism which depends 
critically on the metallicity. It has therefore been suggested that (3 Cep and SPB 
stars should be rare in the Magellanic Clouds which have lower metallicities than the 
solar neighborhood. To test this prediction we have observed the open SMC cluster 
NGC 371 for 12 nights in order to search for (3 Cep and SPB stars. Surprisingly, we 
find 29 short-period B-type variables in the upper part of the main sequence, many 
of which are probably SPB stars. This result indicates that pulsation is still driven 
by the k mechanism even in low metallicity environments. All the identified variables 
have periods longer than the fundamental radial period which means that they cannot 
be (3 Cep stars. Within an amplitude detection limit of 5 mmag no stars in the top 
of the HR-diagram show variability with periods shorter than the fundamental radial 
period. So if (3 Cep stars are present in the cluster they oscillate with amplitudes below 
5 mmag, which is significantly lower than the mean amplitude of (3 Cep stars in the 
Galaxy. We see evidence that multimode pulsation is more common in the upper part 
of the main sequence than in the lower. We have also identified 5 eclipsing binaries 
and 3 periodic pulsating Be stars in the cluster field. 
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1 INTRODUCTION 

f3 Cep and slowly pulsating B stars (SPB) pulsate due 
to the K-mechan ism activated by the metal opacity bump 
|Cox et al.l 1 19921 ). This give rise to p-mode pulsation in (3 
Cep stars and g-mode pulsation in SPB stars. However the 
theoretical standard models do not predict the presence of 
pulsation in (3 Cep and SPB stars in low-metallicity envi- 
ronments such as the Magellanic Clouds. 

Analysis of OGLE data have indeed shown that 
(3 Cep and SPB stars exist i n the Magellanic Clouds 
(Pigulski & Kolaczkowskl |2002| ; iKolaczkowski et al. 1 12004 
2006) though it is not clear if the rate of (3 Cep and SPB 
stars is lower than or equal to the rate in the Galaxy. We 
have therefore observed the open SMC cluster NGC 371 for 
12 nights in order to search for (3 Cep and SPB stars. Ob- 
serving a single cluster instead of the the entire SMC has 
a number of advantages, - e.g. the age and metallicity of 
the cluster can be constrained from fitting isochrones to the 
HR-diagram. 

The metallicity of the SMC has been measured to be- 
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tween Z = 0.001 and Z = 0. 004 fsee iMaeder et al1|l999l . 
and references therein) , while iMiglio et al.l ( 2007al lbf) have 
calculated instability domains of /3 Cep and SPB stars us- 
ing improved opacities and metal abundances which do not 
predict pulsation in (3 Cep and SPB stars for Z < 0.005. We 
should therefore not expect to find any j3 Cep or SPB stars 
in the SMC. 

In order to test this prediction we have chosen to ob- 
serve the cluster NGC 371 as this c luster has been observed 
before by [Kjeldsen fc Baadel (|l994l ) who saw some signs of 
variability, though no clear evidence of pulsation was seen. 

The observations presented in this paper have been 
made with a single telescope. This means that clear mode 
and frequency determination can not be obtained and the 
frequencies can not be used for modeling. Instead this work 
presents a number of candidate SPB stars for followup ob- 
servations with at least dual-side photometry and with spec- 
troscopy. 

The paper is arranged as follows. Section 2 presents 
the cluster. Section 3 describes how we obtained, reduced 
and analyzed the data. The identified eclipsing binaries are 
presented in Section 4, the Be stars in Section 5 and the 
pulsating stars in the upper part of the main sequence in 



2 Karoff et al. 



, . OO \.:v7 vao^.- ry 



. V22 

. -. • .'* ■ I ' 



: V37 '* 



V10- 



...9 « ■ • ® i ? 

^C2] •• . . via . 



V23 



'V35 

' o 



#12 



■:. ■ * .V2$?. G ...-z--. . ■ 
* ^ • ■ -V8 >is ••vfr ' v. - 

; ■ • ■ O 
. •; '■■ vii . .'v. -• : - 



V1 

o 



Figure 1. DFOSC image of NGC371, with the variables marked. 



Section 6. A summary and concluding remarks are found in 
Section 7. 



2 THE TARGET CLUSTER: NGC 371 

NGC 371 (a 2 ooo, 8 20 oo = ll h 03 I °25!0 -72°04 / 40 , /0) is a 
young open cluster in the SMC. IWisniewski fc Biorkmar] 
(2006) estimated a log(age) of 6.7 based on isochrone fitting 
to OGLE data assuming a metallicity of Z = 0.002. As there 
are no high-resolution spectra available for NGC 371 it has 
not been possible to estimate the metallicity autonomously. 

The field of the cluster is shown in Fig. 1 with the iden- 
tified variable stars marked. 
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Figure 2. Light curves of one of the bright stars in NGC 371 
showing the precision level and data sampling. 



BIAS and the flat-fields were indeed stable over the length 
of the observing run. 

The photometric reductions were done using the soft- 
ware package MOMF l|Kieldsen fc Frandseniri992[ ). MOMF 
applies a very robust algorithm combining PSF and aperture 
photometry in semi-crowded fields. 

Some of the light curves showed residuals of systematic 
trends caused by, e.g., air mass, cloud cover, and tempera- 
ture variations (these residuals are sometimes referred to as 
"r ed noise"). We there fore applied the algorithm introduced 
bv lTamuz et all i|2005l ) in order to correct for systematic ef- 
fects. This clearly improved the rms noise level, in particular 
in the brightest stars. 

The resulting light curves had rms noise levels over the 
whole observing run ranging from a few mmag in the bright 
end and up to 50 mmag in the faint end. Some of the light 
curves showed night-to-night drifts, which means that the 
rms noise level in the individual nights was lower than over 
the whole observing run. Fig. 2 shows the precision in a 
typical bright star. 



3 THE DATA AND THE DATA REDUCTION 

The observations were obtained with the DFOSC instru- 
ment at the Danish 1.54-m telescope at ESO, La Silla during 
12 nights in August and September 2005. 763 frames were 
collected in I and 623 in B. The same pointing was main- 
tained during the observations, i.e., the stars were kept at 
fixed positions on the CCD (within a few pixels) during the 
observing run. Because of the crowding in the field the ob- 
servations were always made in focus and the exposure times 
were then adjusted according to the seeing so that only the 
brightest 5 % of the stars were saturated. This resulted in 
exposure times of approximately 20 s in I and 60 s in B. In 
Fig. 2 all the data are shown for one of the bright stars in 
NGC 371, illustrating the time distribution of the data set. 

The CCD images were calibrated using standard proce- 
dures, i.e. a master BIAS was subtracted from each frame 
before the frames were divided by a master sky-flat in each 
filter. The master BIAS was obtained from a large num- 
ber of frames from the whole observing run and the master 
sky-flats were obtained as the median of a large number of 
evening and morning sky-flats. We checked that both the 



3.1 Analysis of differential light curves 

We obtained differential light curves of 6842 stars in the 13.5' 
x 13.5' FOV and we searched for variability in all these light 
curves. The light curves were kept in the instrumental sys- 
tem, but calib r ated UBVI photometry was obtained from 
IZaritskv et~ai1 (|2002h . 

We used two complementary algorithms to search 
for variability; the analysis of variance periodograms 



l)Schwarzenberg-Czernvl 19961) and sim ultaneous iterative 
sine wave fitt ing |Frandsen et al.|[l995h based on the Lomb 
periodogram l|Lomblll976l )7 We searched for variability both 
with and without using statistical weights (Frandscn et al. 
1995). The statistical weights were assigned as the standard 
deviation of all data points separated by less than 15 min- 
utes from a given data point. For the analysis of variance 
approach we made visual examination of all light curves 
which had fit qualities better than 0.9 and for the Lomb 
periodogram approach we made visual examination of all 
light curves where the same peak in the periodogram was 
present in both B and / with a S/N higher than 4 in ampli- 
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Figure 3. Phase light curves of the 5 eclipsing binaries identified 
in NGC 371 



Table 1. Stellar parameters for 5 eclipsing binaries in NGC 371. 
Periods are in days 



ID 



02000 



B - I 



VI 


01 03 04.1 


-72 08 04 


18.74 


0.07 


1.348 


V2 


01 03 34.9 


-72 08 60 


18.64 


2.68 


1.045 


V3 


01 03 40.5 


-72 01 26 


18.86 


-0.17 


0.556 


VI 


01 30 40.0 


-72 02 33 


19.14 


-0.22 


0.463 


V5 


01 03 43.5 


-72 05 11 


18.64 


-0.23 


0.609 



Table 2. Stellar parameters for 3 Be st ars in NGC 371. The id's 
in the second column (ID II) are from IWisniewski fc Biorkmanl 
J2006h . 



ID 


ID II 


«2000 


^2000 


B 


B - I 


V6 


WBBe 43 


01 03 34.5 


-72 06 42 


16.65 


0.35 


V7 


WBBe 110 


01 03 30.5 


-72 03 46 


17.99 


-0.34 


V8 


WBBe 18 


01 03 52.6 


-72 05 39 


15.73 


0.21 



4 ECLIPSING BINARIES 

We have detected 5 eclipsing binaries (phased light curves 
are shown in Fig. 3). VI, V2 and V3 are located on the 
edge of the cluster, which means that they might not be 
cluster members. V4 and V5 are located safely inside the 
cluster and these stars could be important in determining 
the precise distance and age of the cluster. This information 
will again be important when trying to model the excitation 
mechanisms of /3 Cep and SPB stars based on data from this 
cluster. 

V2 also shows some signs of pulsation with a frequency 
of 6.2 c/d, which suggest that this could be a binary sys- 
tem with one of the members being a pulsating star. But 
more photometry will be needed in order to evaluate the 
phenomenon properly, - i.e removing the eclipse from the 
light curve before analyzing the pulsation. 



tude ijBreeer et al.lll993l ). In total we ended up with a list 
of a little over a hundred stars that were selected for further 
analysis. 

In the visual examination, stars were rejected mainly 
because of one of the following three causes: stars were 
placed outside the cluster region, stars were placed close 
to a hot pixel or a bad column, or the variability was not 
the same in the two filters. This gave us the list of the 37 
stars presented in this paper. Light curves of all the variable 
stars will be added to CDS. 

It is not possible to determine if a star is a true clus- 
ter member or not, as we only have two dimensional images 
available. To determine if a star is likely to be a true clus- 
ter member in three dimension a spectroscopic analysis is 
needed. Therefore we can not give a reliable estimate of the 
number of cluster members and some of the identified vari- 
able stars could turn out not to be clusters members. 



5 BE STARS 

Three of the detected var iable stars match with t h e star s 
identified as Be stars by IWisniewski fc Biorkmanl (|2006h . 
They have identified 118 Be stars and 11 candidate Be stars 
in NGC371 based on two-color diagrams of B, V, R and Ha 
photometry. Though the nature of the variability of Be stars 
is believed to be transient, Be stars might also exhibit g- or 
p-mode pulsation. It is therefore clear that these three stars 
can be used to gain understanding of the relation between 
g- and p-mode pulsation and the variability of Be stars in 
low-metallically environments. The amplitude spectra of the 
three Be stars are shown in Fig. 4. The light curves of the 
three Be stars all show coherent variability, which indicate 
that the variability originates from pulsation and not from 
transient events like activity. 
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Figure 4. Amplitude spectra for 3 Be stars in NGC371. Note 
that the scale on the y-axis is different for the three stars. 



Figure 5. Amplitude spectra for three pulsating stars in the 
upper part of the main sequence of NGC 371 



6 PULSATING STARS IN THE UPPER PART 
OF THE MAIN SEQUENCE 

The frequency analysis of the varia ble stars in th e uppe r 
part of the main sequence followed lArentoft et all (2007). 
This means that we made a simultaneous least-squares fit 
to all the frequencies with a S/N higher than 4 in the am- 
plitude spectra of both the I and the B filter data. For 
each frequency we manually inspected the amplitude spec- 
tra, in order to find the frequency value that best described 
the variation in both filters. This step could not be made 
completely objective as daily aliases were present in the am- 
plitude spectra because the observations were made from a 
single site. Therefore some of the detected frequencies might 
be shifted 1-2 c/d from the true oscillation frequencies. The 
uncertainties on the amplitudes and phases were calculated 
based on the error matrix in the least-squares fitting proce- 
dure, though we also used Monte Carlo simulations to check 
the consistency of t hese uncertainties [th is was done with 
the use of Period04 l|Lenz fc Breed l2005h ]. This procedure 
for calculating the uncertainties assumes that we know the 
frequencies, - i.e. that we know which peak reflect the oscil- 
lation mode. The calculated uncertainties does therefore not 
reflect the uncertainty originating from mismatch between 
the true oscillation mode and the daily aliases. 

All the detected frequencies for the pulsating stars are 
given in Table 3 together with color, amplitude and phase 
information, three examples of amplitude spectra of pulsat- 
ing stars are shown in Fig. 5 and phased light curves are 
shown in Fig. 6. 
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Figure 6. Phased light curves for three pulsating stars in the 
upper part of the main sequence of NGC 371 
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Table 3. Stellar parameters for the pulsating stars in NGC 371. Frequencies (y) are in c/d, amplitudes (Ag) are in mmag, phase 
differences (<j>j — <j>g) are in radians and pulsation constants Q are in d _1 . The quoted errors on the amplitude ratios and phase 
differences are based on the error matrix in the least-squares fitting procedure. 
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6.1 Color-magnitude-diagram of NGC 371 



We have used the photometry from lZaritskv et alJ (|2002l ) to 
make a B versus B — I color-magnitude-diagram of NGC 
371 as shown in Fig. 7. The center of the cluster was placed 
on 02000, fcooo = ll h 03 m 25':0,-72°04'40'.'0 and the size was 
selected as a circle with a radius of 5'. The color-magnitude- 
diagram shows the upper part of the main sequence, which 
is almost vertical, and the giant branch. All the identified 
variable stars which are not eclipsing binaries or Be stars 
are located on the upper part of the main sequence, which 
is also where the instability strips for /3 Cep and SPB stars 
would be located if the metallicity of the cluster had been 
larger. Instability analysis of models of hot main sequence 
stars of solar metallicity show a group of mostly p-mode 
pulsators where the longest period is the fundamental radial 
mode and a second group of cooler stars with mostly g- 
modes and longer periods. The first group can be identified 
with the classical j3 Cep stars (early B-type stars) while the 
second group can be identified with the SPB stars (mid B- 



type stars). Models based on revised solar metal abundances 
and different opacities show that the hot boundary of the 
SPB instability strip ov erlaps with the (5 Cep instability 
strip (jMiglio et al.ll2007al lbl). At the metallicity of NGC 371, 
however, no pulsations are excited in the models as already 
mentioned. Therefore theory cannot be used as a guide in 
classification. We will come back to this in the next section. 

We have us e d the Padova evolutionary code described 
m iMarigo et ail (|2007l ) to calculate an isochrone for the 
color-magnitude-diagram of the cluster. The model used so- 
lar heavy elements mixture. We adopted a metalicity of 
Z=0.0 02 and a log(age) of 6.7 from lWisniewski fc Biorkmanl 
J2006I). The distance modulus was chosen to 18.7 from 
ICrowl et all (l200lh and the reddening to E(B - V) = 0.085 
i|Larsen et al.ll2000l ). 



The isochrone is plotted on top of the color-magnitude- 
diagram in Fig. 7. Here it is seen that the isochrone generally 
fits the data nicely. 
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Figure 7. Color-magnitude-diagram of NGC 371 with the iden- 
tified SPB candidates marked. The solid line shows a isochronc 
fit to the observations. 



6.2 Period distribution in NGC 371 



Figure 8. Magnitude-period relation for the candidate SPB 
stars stars together with the fundamental radial period from the 
isochrone. The mode with the longest period in each star is plot- 
ted with large crosses, while the secondary modes with shorter 
periods is plotted with small diamonds. 



In the Galaxy, f3 Cephei stars have periods which are con- 
sistent with mostly p-mode pulsations. The longest period 
is that of the fundamental ra dial mode for which the p ulsa- 
tion constant Q » 0.033 d" 1 |Stankov fc Handle rlliooa ) . To 
classify stars in NGC 371 we need an estimate of Q for each 
star: 

logQ = logPT l/21og(M/M ) - 3/41og(L/L Q ) 

+31og(T/T ). 

where P is the period in days, M the mass, L the luminos- 
ity and T the effective temperature. Since we do not have 
data to estimate the masses, luminosities and effective tem- 
peratures for individual stars, we use these values of the 
isochrone shown in Fig. 7. The result is shown in Fig. 8 
together with all the identified periods as a function of B 
magnitude. There the mode with the longest period in each 
star is plotted with large crosses, while the secondary modes 
with shorter periods are plotted with small diamonds. 

In Fig. 8 it is seen that all stars have periods that are 
longer than the fundamental radial periods. We therefore 
identify all the stars as candidate SPB stars. Classification 
of the stars as bona fide SPB stars would require more time 
series observations to ensure that the periods are correct and 
exclude other sources of variability, such as close binaries. 

If there are (3 Cep stars in the cluster then we can 
assume two things about them; That they are generally 
brighter than the SPB stars and that they have periods 
shorter than the fundamental radial period. We do not de- 
tect any pulsating stars in this domain than fulfill the de- 
tection criteria given in Section 3.1 and therefore we do not 
detect any candidate (3 Cep stars in NGC 371. 

In order to constrain an upper detection limit of (3 Cep 
stars in the cluster we measured the highest peak in the 
amplitude spectrum between 8 and 20 c/d for all the stars 
brighter than 15.8 in B. A frequency range from 8 to 20 
c/d is in agreement with the frequency range in which the 
f3 Cep stars are expected for this cluster (see Fig. 8). For 
the brightest stars we could also have included periods up 
to 0.4 day, but this gave problems with the 1/f noise in the 
amplitude spectra of some of the stars. 

In Fig. 9 we have plotted the amplitude of the high- 



est peaks as a function of magnitude together with the 
mean V oscillation ampli tude of galactic j3 Cep stars from 
IStankov fc Handler! (2005). It is seen that our detection limit 
is significantly lower than the mean V oscillation amplitude 
of galactic /3 Cep. If (3 Cep stars are present in the cluster 
they have oscillation amplitudes significantly lower than the 
galactic j3 Cep stars. 

Though we do not detect any candidate [3 Cep stars 
some of the stars do show variability between 8 and 20 c/d, 
but not with a S/N higher than 4 in both filters. In Fig. 9 we 
have marked the peaks which have S/N higher than 3.5 in 
B and in Fig. 10 we have plotted the amplitude spectrum of 
one of them. The detection criteria of S/N higher than 4 in 
both filters is generally conservative and for longer regularly 
sampled data sets it is too high as the observations in the 
two filters are independent. But for the present data set it 
seems as a good indication of pulsation. More photometric 
data is needed in order to classify these stars as candidate ft 
Cep stars. On the other hand if these stars are /3 Cep stars 
they would be very interesting targets as they would have 
shorter periods and lower oscillation amplitudes compared 
to the galactic j3 Cep stars. 

Though periods below 0.1 day are not common in (3 Cep 
stars they have been seen in the very young cluster NGC 
6231 (|Balona fc Shobbrook|[T983l ; lArentoft et al.ll200ll ). 

6.3 Amplitude ratio and phase difference 

We have also calculated and plotted the amplitude ratios 
and phase difference between the B and I filters (Fig. 11) in 
order to identify any systematic effe cts that could be com - 
pared with theoretical models as in [ Arcntoft et alj (I2007T) 
or used for mode identification (|Hevnderickx et Uhl Il994h . 
Though we do not find any systematic effects in phase 
difference versus amplitude ratio within the uncertainties 
(Fig. 11) we do see that the main part of the observed fre- 
quencies cluster around an amplitude ratio between 0.5 and 
1.0 which is expected fo r oscillations in f3 Cep s tars with 
degrees between and 2 (jHevnderickx et alii 19941 ). 

The I/B amplitude ratio can also be used to exclude 
other causes of the variability than pulsation (e.g. ellipsoidal 
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Figure 9. Amplitudes of the highest peak as a function of mag- 
nitude for the stars brighter than the identified candidate SPB 
stars. The diamonds mark the peaks with a S/N higher than 3.5. 
The horizontal line at 17 mmag shows the mean oscillation am- 
plitude of /3 Cep stars in the Galaxy iaV . The arrow marks the 
star plotted in Fig. 10. 

variables) as t he I/B amplitude r atios are general below 
unity (see, e.g., lArentoft et "al]|2007| , and references therein). 

6.4 Multimode pulsation 

The last thing to note about the frequencies is that all the 
multimode pulsators are found among the pulsating stars 
in the upper part of the main sequence (except for one). 
One reason for this could be that these stars are brighter, 
which in general gives a higher S/N, but some of the fre- 
quencies that we detect in the stars in the lower part of the 
main sequence have such a high amplitude that multimode 
pulsation should have been discovered if it was present (see 
Table 3). 



7 SUMMARY AND CONCLUSIONS 

We have identified 5 eclipsing binaries, 3 periodic Be stars 
and 29 candidate SPB stars in NGC371. 

The results indicates that excitation of oscillations in 
SPB stars is more common in low metallicity environments 
such as the SMC than predicted by standard stellar mod- 
els. As a possible extension to the standard stellar models 
to account for this discrepancy could be to include local 
iron enhancement b y diffusion and radiative accelerations 
l|Miglio et al.ll2007d\ 

If there are /3 Cep stars in the cluster they would gener- 
ally be brighter than the SPB stars and have periods shorter 
than the fundamental radial period and we do not detect 
pulsation in any stars in this domain with an upper limit on 
the oscillation amplitudes of 5 mmag. 

Though it is not possible to calculate a reliable num- 
ber for the fraction of SPB stars in NGC 371 the absolute 
number of pulsating stars seems to b e high compared to 
the fr action obtained for the Galaxy bv lStankov fe Handler! 
(120051) . 

We see evidence that multimode pulsation is more com- 
mon in upper part of the main sequence than in the lower. 
We have also identified periodic pulsation in 3 stars that 
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Figure 10. Amplitude spectrum in B of one of the stars showing 
signs of /3 Cep variability between 8 and 20 c/d with a S/N below 
4. 

were previously identified as Be stars. These stars can used 
to study the relation between the Be phenomenon and j3 
Cep and SPB stars in low-metallicity environments. 

The results presented here strongly confirm and increase 
our interest in NGC 371 as an excellent laboratory for j3 Cep 
and SPB stars. 
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